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The shear viscosity of a two-component Fermi gas in the normal phase is calculated as a function 
of temperature in the unitarity limit, taking into account strong-coupling effects that give rise to 
a pseudogap in the spectral density for single-particle excitations. The results indicate that recent 
measurements of the damping of collective modes in trapped atomic clouds can be understood in 
terms of hydrodynamics, with a decay rate given by the viscosity integrated over an effective volume 
of the cloud. 
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I. INTRODUCTION 

Strongly interacting Fermi systems play a central role 
in physics over a vast range of energies, from cold atoms 
over condensed matter systems to quark-gluon plasmas. 
For atomic gases, the regime of strong interaction is 
reached by the use of Feshbach resonances at which the 
scattering length diverges. This so-called unitarity limit 
has been studied experimentally through the expansion 
of a two-component Fermi gas [l[ and by measuring 
its collective modes 0, These experiments indicate 
that under certain conditions the dynamic properties of 
atomic gases in the unitarity limit are well described 
by hydrodynamics, both in the superfluid and in the 
normal phase. Related results have been reported for 
the strongly interacting quark-gluon plasma produced in 
heavy- ion collisions at RHIC [j, Q ■ 

In this paper we shall carry out a quantitative analysis 
of the hydrodynamic damping for the normal phase in 
the unitarity limit and compare our results to the mea- 
sured rate of decay of the collective modes @, Q . In the 
hydrodynamic limit, the rate of decay is related to the 
shear viscosity integrated over the volume of the trapped 
atomic cloud. As we shall see, under the given experi- 
mental conditions, hydrodynamics applies only in a lim- 
ited temperature region above the superfluid transition 
temperature T c . Furthermore, hydrodynamics necessar- 
ily fails in the outer parts of the atomic cloud, where 
the density is low and the mean free path therefore long. 
Since for a classical gas the viscosity is independent of 
density, one must introduce an explicit cut-off in the spa- 
tial integration of the viscosity, as shown in Ref. @. 

In the unitarity limit, we can use a dimensional argu- 
ment [7[ to write 77 as 

7? = nha(T/T F ). (1) 

Here a is a dimcnsionless quantity which can only de- 
pend on temperature through T/Tp. The Fermi temper- 
ature is Tp — kp/2m, where kp is the magnitude of the 
Fermi momentum given by k F = (37r 2 n) 2 / 3 with n being 
the density of the gas (with the exception of (P), we use 
units such that h = ks = !)■ Our aim is to obtain an ap- 



proximate expression for the universal function a(T/Tp) 
which will allow us to compare theory with experiment. 

The present work is a continuation of two previous pa- 
pers H Q , in the following referred to as I and II respec- 
tively, on the damping of collective modes in Fermi gases. 
Before we describe the results of our present calculation 
we shall therefore summarize the approach taken in I and 
II, and indicate their limits of applicability. 

In I we employed a Boltzmann equation for the fermion 
distribution function /(r,p,t), taking into account the 
dependence of the scattering cross section on the energy 
in the relative motion of two particles, as well as the effect 
of the mean field in the streaming terms of the Boltz- 
mann equation. The collective mode frequencies were 
calculated by taking moments of the Boltzmann equa- 
tion, which introduced the spatially-averaged viscous re- 
laxation rate as the effective collision rate entering the 
imaginary part of the (complex) mode frequencies. At 
low temperatures, well below the Fermi temperature, the 
validity of this approach is restricted to the limit of weak 
coupling, kp \a\ <C 1, where a is the scattering length and 
kp is the magnitude of the Fermi wave vector. In this 
limit the method is accurate within a few percent, since 
the viscous relaxation rate used in I is obtained from a 
trial function which yields viscosities that differ by only 
a few percent from those obtained from exact solutions 
to the Boltzmann equation at both low and high tem- 
peratures. At temperatures well above the Fermi tem- 
perature the Boltzmann approach used in I is accurate 
for any value of a, including the unitarity limit in which 
the cross section is proportional to the inverse of the en- 
ergy in the relative motion. We shall demonstrate this in 
detail in App. |A"1 

When kp\a\ becomes comparable to or larger than 
unity, one enters the regime of strong coupling, in which 
perturbative approximation schemes such as the Boltz- 
mann approach can no longer be trusted at tempera- 
tures comparable to or less than the Fermi temperature. 
Progress in understanding the viscosity of such strongly 
coupled Fermi systems must necessarily rely on an in- 
terplay between experiment and theory, since there is 
no small parameter available for a perturbation expan- 
sion that could yield firm theoretical predictions. For an 
atomic gas close to a Feshbach resonance we explored in 
II the influence of the medium on the scattering cross 
section, which in I was taken to be its value in vacuum. 
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Due to Fermi blocking of the pair states into which the 
molecular state can decay, the lifetime of the resonant 
state was found to be significantly increased, leading to 
a corresponding increase in scattering rate (and hence a 
decrease in viscosity) close to the superfluid transition 
temperature T c . For a uniform gas the calculated viscos- 
ity just above T c was found to be reduced by the factor 
7.5 compared to the value obtained in I by use of the vac- 
uum scattering matrix. For a trapped gas the difference 
was less pronounced: the thermal relaxation rate, which 
is closely related to the inverse viscosity, was found near 
T c to be 3.6 times the value obtained using the vacuum 
scattering matrix. 

The assumption underlying the approach taken in II 
was that the main effects of the interaction in the strong- 
coupling limit arose through a modification of the scatter- 
ing cross section, while strong-coupling effects that lead 
to spectral broadening of single-particle excitations were 
not taken into account. Put in different terms, only the 
collision term in the kinetic equation was modified by 
taking into account the medium effects mentioned above, 
while the streaming terms were assumed to be unaffected 
by interactions. In the present paper we abandon this as- 
sumption and consider specifically the role of the pseudo- 
gap occurring in the spectral function of single-particle 
excitations |l0(]. The presence of the pseudogap in the 
normal phase influences the Bragg spectrum observed 
when an atom absorbs a photon from one laser beam and 
emits a photon into another, resultin g in a change of the 
energy and momentum of the atom [111 ] . In the normal 
phase, at the unitarity limit, the pseudogap was found 
to cause a significant suppression of the low-frequency 
Bragg spectrum. 

The use of a Boltzmann equation as in I and II im- 
plicitly assumes the existence of quasiparticles with a 
definite energy-momentum relationship. When the spec- 
tral functions broaden, the quasiparticles are less well- 
defined, and it therefore becomes relevant to investigate 
the effect of this broadening on the transport properties 
of the gas. Ideally, one should derive a kinetic equation 
that takes all strong-coupling effects systematically into 
account, but due to the lack of a small parameter in the 
strong-coupling limit this would be far too ambitious an 
undertaking. Our aim here is more modest: to com- 
pare results for the viscosity in the presence and absence 
of spectral broadening in order to gain insight into its 
quantitative significance, and to use this together with 
the results of I and II to construct an approximate for- 
mula for a that allows for an explicit comparison with 
experiment. 

A main result of the paper is the calculation in Sec. 
HTl of the shear viscosity from a Kubo formula, which al- 
lows one to take into account the presence of the pseu- 
dogap in the spectral density of states for single-particle 
excitations. Since our results are based on a ladder ap- 
proximation to the self-energies, they cannot be expected 
to be quantitatively accurate, but as we shall see, our 
present results are quite close to those of paper II at low 



temperatures. At high temperatures, however, the Kubo 
expression gives results that lie significantly below the 
classical result obtained from the Boltzmann equation. 
Since the latter is essentially exact for all values of a in 
the classical limit, as we shall demonstrate in App. [A] 
we construct an approximate formula for a which inter- 
polates between the low-temperature result and the ex- 
act high-temperature limit. This interpolation formula 
is then used for comparison with experiment in Sec. IIII1 
where the decay rate is related to the viscosity integrated 
over an effective volume of the trapped gas. The resulting 
agreement with experiments 0, Q that were carried out 
at two very different frequencies indicates that the inter- 
polation formula provides a qualitatively correct picture 
of the physics involved in the viscosity of a strongly in- 
teracting Fermi gas. The calculations also illustrate how 
information on the viscosity of strongly interacting Fermi 
gases can be extracted from measurements of the damp- 
ing of collective modes. 



II. THE SHEAR VISCOSITY 

Consider the shear viscosity 77 of a uniform, two- 
component Fermi gas of atoms with mass to in the normal 
phase. There is no interaction between atoms in the same 
internal state whereas the interaction between atoms in 
the two different internal states is characterized by the 
s-wave scattering length a. Unitarity means that we take 
kp\a\ — ► 00. The shear viscosity relates the momentum 
current density to the gradient in flow velocity u x (y) 
according to n^j, = —rjdu x /dy. 

The Landau-Boltzmann approach assumes well- 
defined quasiparticle excitations. However, with strong 
interactions present the spectral functions may become 
significantly broadened and the quasiparticles therefore 
ill-defined. Close to T c in the normal phase, the spectral 
weight is found to be suppressed near the Fermi surface 
resul ting in a double-peak structure of the spectral func- 
tion [lO , [ll I . This suppression is referred to as the pres- 
ence of a pseudogap. In order to investigate the impor- 
tance of the pseudogap we turn to the Kubo formalism, 
which allows for a treatment of these effects. 

The velocity field u x {y) gives rise to a perturbation 



H' = -to / d 3 ru x (y)j x (r) = 



1 du x 
iuj dy 



d 3 rfl xy (2) 



with j being the current density operator and tl xy the 
momentum-current density operator. To obtain the sec- 
ond equality in ([2]), we have used momentum conserva- 
tion mdji/dt = —dUij /drj (i and j = x,y,z) taking 
j(r,t) — exp(— iwi)j(r). We can now calculate H xy in- 
duced by H' within linear response. The shear viscosity 
is then obtained by taking the limit lu — > 0: 



T] = — lim Ima(u)/w 



(3) 
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equal to their equilibrium values, 



FIG. 1: 3 and the atomic propagator G in the ladder approx- 
imation for a broad resonance. E is the self energy and V is 
the scattering matrix. 



with 



3(w) = -* / d a rdte iut e(t)([tL xy (r,t),IL xy (0,0)]). (4) 



The approximation for 3 shown in Fig. [T] yields 



V = 



15m 2 



d 3 p 4 
P 



(27T) 



2n oe 



(5) 



Here, A(p, e) = —2lmGn(p, e) is the spectral function for 
the atoms with Gr being the retarded Green function 
and /(e) = [exp(e/T) + A relativistic version of ([5]) 
has been used to calculate r\ for quark-gluon plasmas us- 
ing a phenomenological ansatz for A(p,e) Here, we 
use a microscopic theory for A(j>, e) as explained below. 

To obtain A{p, e), we calculate the Green function us- 
ing a multichannel BEC-BCS crossover theory based on 
the ladder approximation for the thermodynamic poten- 
tial and the atom self energies. This is the minimal ap- 
proach which includes the correct two-body physics lead- 
ing to the Feshbach resonance [HI, 0, EH • The structure 
of the theory is shown in Fig. [1] and is described in detail 
in Refs. 
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Il5j . We take a broad resonance for which 
the multichannel theory becomes equivalent to the orig- 
inal single channel BEC-BCS crossover theory [l6| for 
most observables. The spectral functions are found nu- 
merically to obey the sum rule / A(p, oj)dui — 2tt to a 
very good approximation. It should be noted, however, 
that the approximation leading to © is not conserving. 
To obtain a conserving approximation, we need to solve 
an integral equation for S jl7j ]. However, the present 
analysis is already heavy numerically since one needs to 
integrate over the atom self energies which themselves 
are evaluated numerically; a conserving approximation is 
thus beyond the scope of the present paper. 

It is instructive to compare the Kubo approach with 
the relaxation-time approximation to the Boltzmann 
equation, since these yield identical results at high tem- 
peratures. In the relaxation-time approximation, the col- 
lision integral of the Boltzmann equation given in Eq. 
(13) of Ref. § becomes /[/] = Sf/r, where Sf is the 
deviation of the distribution function from local equilib- 
rium. The relaxation time r(p) is obtained by setting the 
distribution functions for particles in states other than p 



1 



d 3 Pl 



dn 



da |p - pi| 



m 



r{p) J (2tt) 3 J ~~~dtt 

1 - /(&) 



(6) 



Here da/dQ. is the differential cross section for the scat- 
tering of two particles with incoming momenta p and 
Pi to outgoing momenta p' and p^. The correspond- 
ing energies are £ p = p 2 /2m — fi etc., and Q is the 
solid angle of (p' — p[)/2 with respect to (p — Pi)/2. 
It is straightforward to find Sf from the Boltzmann 
equation and thus the momentum current density from 
U xy — 2(2-7r) -3 J d 3 p(p x p y /m)Sf, where the factor of 2 
is from the two internal states. We obtain for the shear 
viscosity 

d 3 p *df 



15m 2 



W at, 



(7) 



The factor 1/15 arises from the angular average of p 2 p 2 . 

The Kubo formula © reduces to the relaxation- 
time approximation © in the limit of weak interaction 
with narrow spectral function peaks. To see this, put 
A(p,e) = 2r p /[(e - £ p ) 2 + T 2 } with T p = -ImE(p,&) 
the imaginary part of the atom self energy. Using 
A 2 -> 27r<5(e - £p)/ r p for r p -> 0, © reduces to © with 



-Hp) 



-2ImS(p, ^ p ). It can furthermore be shown 



that the ladder approximation for lmS(p, £ p ) yields (J6J) 
when medium effects are ignored [l5l |. The Kubo for- 
mula for the viscosity ([5]) thus reduces to ||7J) in the limit 
of weak interactions. 

In Fig. [21 we plot a(T/T F ) in for a gas in the uni- 
tarity limit with T > T c . For the numerical calculations, 
we have chosen parameters corresponding to a resonant 
interaction with kp\a\ = 25 3> 1 and a negligible ef- 
fective range. The critical temperature is found from 
the divergence of the scattering matrix at zero total mo- 
mentum (the Thouless criterion) yielding T c « 0.26Tp 
in good agreement with other BEC-BCS crossover re- 
sults [H, [HI]. The Kubo result © approaches © for 
T ^S> Tp. This is to be expected since medium effects arc 
negligible in the classical limit. For T ^> Tp, 17 cx T 3 / 2 
which can be seen in the relaxation-time approximation 
from ©-© which give 77 = nTr with r cx T 1 / 2 , re- 
sulting in a = 1.1(T/T F ) 3 / 2 . At low T, the difference 
between © and © is significant; for T = T c , © yields 
a = 0.4 whereas © predicts a — 0.2. This reduction is 
due to strong-coupling medium effects. The imaginary 
part ImS is large leading to a significant damping of the 
quasiparticles. Close to T c a pseudogap opens up and 
the effective density of states at the Fermi surface is sup- 
pressed, leading to a reduction of the viscosity. 

We also plot 77 obtained from a variational solution 
to the Boltzmann equation [1, [§], which yields a lower 
bound on the viscosity, given by 



45tt 3 / 2 / T\ 
64^2 \Tf~J 



3/2 



T 

= 2.77 I — 



3/2 



(8) 



4 




FIG. 2: (Color online) The viscosity in units of nh for T >T C . 
The dashed-dotted line is the Kubo formula (JjsJ, the dashed 
line is the Boltzmann equation result in the relaxation-time 
approximation {2} , the solid line is the variational solution of 
the Boltzmann equation with the medium cross section [jj, 
and the dotted line is the high-T result (JSj) . 



This lower bound deviates by less than two percent from 
the exact result, as we shall demonstrate in the Ap- 
pendix, where we investigate the leading correction to 
the lowest-order variational result for a general form of 
the cross section. In the unitarity limit the leading cor- 
rection is found to increase the lower bound (|SJ) by 3/190 
or 1.58 %. For comparison we consider in the Appendix 
also the weak-coupling limit where the cross section is a 
constant, independent of the relative momentum of the 
colliding particles. In this case the leading correction is 
found to increase the lower bound by 3/202 or 1.49 % in 
precise agreement with the known result for hard-sphere 
molecules [19j . 

Since medium effects are negligible for T 3> Tp, we 
conclude that Eq. (jHJ is a very accurate expression for 
•q at high temperatures in the unitarity limit. At low T 
however, we saw by comparing ([5]) and j7]) that medium 
effects are important. Compared to the Kubo result, 
the variational solution includes medium effects only in 
the sense that the medium scattering rate r _1 (p) = 
— 2Im£(p, £ p ) is used in the Boltzmann equation implic- 
itly assuming ImE <C fy whereas ReS is neglected 
On the other hand, the variational nature of the solution 
corresponds to approximately summing diagrams beyond 
the approximation in Fig. [1] leading to the Kubo result 
([5]). Comparing the two approximations for r\, we see 
that the variational solution of the Boltzmann equation 
obtained in Ref. [§] agrees reasonably well with the Kubo 
result for low T. Note that since in the unitarity limit 
T c //j,(T c ) « 0.6 is rather large, the T~ 2 divergence in ij 
for T — > due to Fermi blocking is not seen for T > T c 
in Fig. 1 



III. COMPARISON WITH EXPERIMENT 

When conditions are hydrodynamic, the attenuation of 
a collective mode is related to the viscosity. We now use 
our calculated viscosity to interpret the measured 0, Q 
damping of the collective modes in an atomic gas trapped 
in a potential of the form 

V(x, y, z) = ^(lo 2 x x 2 + u; 2 y y 2 + lu 2 z z 2 ). (9) 

The attenuation T of a collective mode is defined in terms 
of the amplitude decay of the density oscillations given 
by one half the rate of change of mechanical energy, i.e. 

p _ (E mcc h)t (10) 
2(-Emcch)t 

where (. . .)* indicates the time average over a period of 
the cycle [6J. The modes we examine are characterized 
by a velocity field v = (ax, by, cz). Following Ref. 0], we 
obtain 

r 2(a 2 + b 2 + c 2 - ab- ac-bc) f d 3 rri(r) | ( 
3m J d 3 r n(r)(a 2 x 2 + b 2 y 2 + c 2 z 2 ) 

for the damping. Here n(r) is the density of atoms 
in the trap. In the unitarity limit, ([I]) gives rj(r) — 
n(v)Tia(T/T F (r)) with T>(r) = [3n 2 n(r)} 2 ^ /2m. We 
parametrize a for T >T C in the form 

a = a() +a 3/2 (T/T F ) 3 / 2 (12) 

where ao = —0.2 and a.3/2 — 2.77 are numbers chosen to 
fit our numerical results for 77 in Sec. Ull 

The ao-term yields a spatial integral of n(r) in (|TT|) 
whereas the a 3 / 2 term gives a spatial integral of a con- 
stant since the viscosity in the classical limit is indepen- 
dent of density. The integration must be cut off near the 
edge of the cloud where the density is small and hence 
the mean free path long, resulting in the breakdown of 
the validity of hydrodynamics. We adopt the procedure 
described in Ref. 6] that the hydrodynamic description 
holds out to a distance given by the condition that an 
atom incident from outside the cloud has a probability 
of no more than 1/e of not suffering a collision. Since 
the density is well approximated by its classical value 
near the edge we obtain as in Q a cutoff distance that 
depends weakly on the cross section a. In the unitarity 
limit this cross section is a = CA^ , where At = \/2-k jmT 
is the thermal de Broglie wavelength and C is a numerical 
constant of order unity. 

In the experiments 0, the trap is very elongated 
with uj z <C uj x ,ujy. The transverse frequencies lu x and 
u) y differ by 10-20%, but to leading order this ellipticity 
can be taken into account for the mode frequencies by 
considering a cylindrically symmetric trap 

V(x,y,z) = ^l(x 2 + y 2 + \ 2 z 2 ) (13) 
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FIG. 3: (Color online) The damping of the transverse mode. 
The x 's are experimental values from [l| and lines are theory. 



with lo± = ^/uj x LOy. To model the experiments, we there- 
fore consider such a trap in the following with A <C 1, 
which results in a separation of the transverse and longi- 
tudinal motion. Hence, a — b and c = for the transverse 
mode, while a = b = and c ^ for the axial mode. 
In the classical limit with r\ given by ([5]) , Eq. yields 



r = i.osx^n- 1 / 3 



U!± 



/(A, r ), (14) 



where /(A, To), given in Ref. [(J, is an angular av- 
erage arising from the integration over volume, with 
to = crn(0)y/T/2m/2Xuj±. The total number of trapped 
atoms is N. For the axial mode, one obtains a similar 
expression for T/ui z with A 2 / 3 replaced by A 5 / 3 , while the 
numerical constant in front is 3.1. 

We shall now compare our calculated damping rates 
with experiments on trapped 6 Li atoms [liH- Measure- 
ments of the damping of the axial mode 2fl yielded the 
value r /uj z 0.0045 in the unitarity limit for a trap with 
A = 0.030 and N = 4 x 10 5 . The temperature, however, 
was not known, and we cannot therefore make a direct 
comparison to our calculated damping. To obtain an esti- 
mate, we use the classical result (| 14[) , which yields values 
in the range 0.004 < T/u z < 0.007 for 0.3 < T/T F < 0.6 

{C = 1) - n 

The Duke experiments [3J| allow for a more direct com- 
parison, since information on the temperature is avail- 
able. In Fig. [31 we plot the observed damping rate of 
the lowest transverse collective mode. The temperature 
is determined by fitting the observed density profile to 
an ideal gas profile with a rescaled Thomas-Fermi ra- 
dius R^p = ^TTPR^f where R 2 TF = 2(3AA) 1 / 3 a 2 L 
and a7 2 = moj\ H- Here, (1 + (5) is the parameter 



= mu>± 3]. Here, (1 
used conventionally in the unitarity limit, equal to the 
ratio between the ground-state energy and that of a non- 
interacting Fermi gas. In the classical regime, the fitted 
temperature T is related to the physical temperature by 
= T/Tp. We have used this relation for all 
T (with (3 = —0.5), since the profile is approximately 
classical above the critical temperature [20l |. 



We plot the calculated damping rate from (fTTj) for two 
different values of C of order 1. The spatial integrals 
are performed with a density profile corresponding to 
an ideal gas with temperature T, Thomas-Fermi length 
Rtf, N = 2 x 10 5 , and A = 0.045. We also plot the classi- 
cal limit (|14p . As expected from Fig. [21 we see that the ef- 
fects of the medium are small except for low T. Note that 
our theoretical curves are not reliable in the region where 
they predict a decreasing damping with increasing tem- 
perature, since this reflects the breakdown of hydrody- 
namic behavior: for larger T the cloud becomes less dense 
and the cross section a also decreases as 1/T, imply- 
ing that the mean free path increases with temperature. 
To estimate the temperature above which the behavior 
ceases to be hydrodynamic, we compare the mean free 
path £(0) ~ l/n(0)<7 in the center to the spatial extent 
of the cloud in the transverse direction. In the classical 
limit £(0) - R x R y R z /Ncr, where R t = (2ttT '/muf) 1 / 2 . 
The condition £(0) < R± implies that T < u^XN) 1 / 2 . 
For the experimental parameters, this means that the 
gas is hydrodynamic for T <C 5/xK~ 2Tp and there is a 
limited temperature range in which we can compare the 
theory to the measured damping rate, since our calcu- 
lations only apply to the normal phase. The fact that 
the observed mode remains hydrodynamic with a small 
damping for T — > indicates superfluidity [2l[ . 

The results for the damping yield the correct order of 
magnitude for two experiments measuring at two very 
different frequencies. This cannot be obtained simply by 
adjusting the parameter C since the results depend only 
weakly on C and we have used C ~ 1 in both cases. 



IV. DISCUSSION AND CONCLUSIONS 

The shear viscosity of a normal Fermi gas in the unitar- 
ity limit was analyzed as a function of temperature. For 
high temperatures where one can perform systematic cal- 
culations, we demonstrated that a variational solution to 
the Boltzmann equation yields a value of rj which is accu- 
rate to better than two percent. At low temperatures, the 
role of the strong-coupling pseudogap effect was analyzed 
by calculating the viscosity within the Kubo formalism. 
We showed that the pseudogap reduces the viscosity sig- 
nificantly since it suppresses the density of states at the 
Fermi level. While we stress that our calculations of the 
viscosity are approximate in nature for T c < T < Tp, it 
is interesting that the Kubo approach yields values that 
are close to those obtained from the Boltzmann equation 
(with a medium cross section) in this regime. This sug- 
gests that the main effects of the medium can be taken 
into account by using a medium cross section in a Boltz- 
mann approach to the transport properties of the atomic 
cloud in the unitarity limit. Based on these high and 
low temperature results, we constructed an interpolation 
formula for the viscosity as a function of temperature at 
unitarity. This formula was used to analyze two experi- 
ments on the decay of collective modes in terms of viscous 
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damping. In performing this analysis, it was crucial to 
introduce an explicit cut-off in the spatial integrations 
since hydrodynamics necessarily fails in the outer parts 
of the cloud, where the density is low. We concluded 
that the hydrodynamic approach of viscous damping ac- 
counts reasonably well for the experimental observations; 
this holds for both the longitudinal and the transverse 
modes for which the observed damping differs by an or- 
der of magnitude. It would be very valuable to be able to 
compare theory and experiment at higher temperatures, 
where the behavior should be collisionless. In this limit 
the mode frequencies should approach twice the oscillator 
frequencies and the damping rate become proportional to 
the inverse relaxation time rather than the viscosity. 

We are grateful to J. E. Thomas for very helpful cor- 
respondence. We also thank T.-L. Ho and C. J. Pethick 
for discussions. 



APPENDIX A: THE HIGH-TEMPERATURE 
LIMIT 

In this appendix we derive the leading correction to 
the lowest-order variational result for the shear viscosity 
of a classical gas. As is well known (see e.g. paper I 
for details) , the properties of the collision integral in the 
Boltzmann equation allow one to derive a lower bound 
on the viscosity of the form 



V 



> 



{U,Xf 

{u.Huy 



(Al) 



where (...,...) denotes a suitably defined scalar product. 
Here X denotes the inhomogeneous term in the Boltz- 
mann equation while U is the trial function and the pos- 
itive, semi-definite integral operator H represents the col- 
lision term. Since we are only interested in determining 
the relative correction to the high-temperature viscosity 
arising from an improved trial function we use units such 
that 2m = ksT = 1 in order to simplify the presentation. 
It is convenient to consider the equivalent upper bound 
on the inverse viscosity, and we shall therefore seek to 
minimize the functional F given by 



F 



(U 1 HU) 
(U,X)* 



(A2) 



by varying the trial function U. The lower bound on the 
viscosity given in ([5]) corresponds to the choice U = X. 

In order to improve this bound we derive the minimum 
value of F for a trial function U given by a variable linear 
combination of the functions U\ and U2, 



U = 7 E7i + c(l - 7)C/ 2 , 



(A3) 



where 7 is a parameter to be varied, while c is a constant. 
We choose c such that (U, X) is independent of 7, 



(Ui,X) 
(U2, X) 



(A4) 



with which we obtain (U, X) — (Ui,X). The numerator 
in (|A2|I is a quadratic form in 7, which is readily mini- 
mized resulting in the value F min . We are interested in 
the relative difference between F m { n and the value F 1= \ 
of F for 7 = 1. Consequently we define 5 according to 



F ■ 

r -, ± mm 
= 1 



(A5) 



where F 1=1 = (C7 X , HU X ) / '{U x , X) 2 : . The integral op- 
erator H is symmetric and therefore (Ui,HU2) — 
(U2,HUi). We define the quantities h\2 and by 



{U 1 ,HU 2 ) 2 (U 2 ,HU 2 ) 
n>\i = c— — and /122 = c 



'{Ui,HU x ) ' {UuHUxY 

in terms of which 8 assumes the form 

s= (i-M 2 

1 + h 2 2 - 2/112 ' 



(A6) 



(A7) 



Note that both h\2 and /122 are independent of any con- 
stant multiplying J7i or U 2 - We take Ui(p) = (p 2 z —p 2 /3) 
(corresponding to U\ oc X) and [^(p) = [p 2 z —p 2 /3)p 2 . 
We also define 

A, = \ (Ui(p) + Ui( Pl ) - U t (p') - [/i(pi)) (A8) 

for i = 1,2. It is convenient to introduce relative and 
center-of-mass momentum variables according to 

Q 



P = q+y Pi 



q+§; P ' = q' + §, Pl = q' + §- 

(A9) 

Since the scattering is elastic, we have q 2 = q' 2 . We 
obtain 



A 2/2 

Ai = q z - q z 



and 



A -j = {ql — q'z) ( 9 2 + -7- ) - -q-Q ('/■/?■- 



q'-Q 



(A10) 

q <r 



q' • Q q' z Q 



(All) 



In order to determine h\2 and /122 we first calculate the 
angular averages (Af), (AiA 2 ) and (A 2 ) by integrating 
over the directions of each of the vectors q, q' and Q. 
We get 



(A?) = 



together with 



(A12) 



(A13) 



and 



.- 0¥ ( 9 ^) +w (2 M (A14) 



7 



The final integrations over the magnitude of the total and 
relative momentum are given by 

/>oo />OC 

/ dqq 2 qa(q) / dQQ 2 e~ 2 ^ ^ ' 2 {. . .), (A15) 
Jo Jo 

where the exponential factors arise from the product of 
the equilibrium distribution functions /(Cp)/(Cpi) an d 
<j(q) is the g-dependent cross section. The additional 
factor of q in the integrand of (|A15|) is due to the relative 
velocity occurring in the collision integral. Putting these 
results together and using that c = r(7/2)/T(9/2) = 2/7 
we finally get 

*»=?(H) (Aw) 

and 

h 22 = ±( I -* + 7 - I J> + ™), (A17) 
22 49 \I 4 2 / 4 48 J ' v ; 



where we have defined the integrals I n by 

/>oo 

In= / dqq n+ Za{q)e~ 2c > 2 . (A18) 
Jo 

In the weak-coupling limit a is independent of q and we 
obtain Ig/14 = 5 and Iq/Ia — 2, resulting in h\ 2 = 15/14 
and h 22 = 8 77/ 588. When inserted into (|AT|) these val- 
ues yield 5 = 3/205 in precise agreement with (l9| . since 
(1 - (5)- 1 = 1 + 3/202. As shown in Ref. [H the conver- 
gence is very fast, the next-order term yielding a further 
correction of only a tenth of a percent. For resonant scat- 
tering a oc q~ 2 , in which case 1%/Ia = 3 and I^/Ii = 3/2, 
resulting in /112 = 13/14 and h 22 = 697/588, which yields 
S = 3/193 or correspondingly (1 — d 7 )^ 1 = 1 + 3/190. 

We can thus safely conclude that the expression ([H]) is 
accurate to better than two percent at all values of a, 
including the unitarity limit a — > 00. 
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